Phase-field study for the sphtting mechanism of coherent misfitting precipitates in 

anisotropic elastic media 
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Coherent misfitting precipitates in elastically stressed media such as 7/ particles in nickel-based 
super-alloys show various splitting patterns such as doublets, quartets, or octets due to their misfit 
strain energy. While it is an interesting instability phenomenon defying conventional surface ther- 
modynamics, its mechanism is not completely clear. Through a phase-field study upon the splitting 
behavior and morphological evolution of coherent precipitates, we show that an interface instability 
driven by elastic anisotropy and a diffusion field can generate elastically induced splitting during 
difTusional phase transition. Particle splitting is triggered by interface grooving which advances by 
penetrating grooves into the interior of the particle. The sequential evolution of shapes during the 
splitting process is in good agreement with previous experiments. 
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Elastic stresses arising during solid state phase tran- 
sitions in some alloys can lead to microstructural fea- 
tures such as precipitate alignment, cuboidal precipitate 
shapes, and particle splitting. A wealth of experimental 
and theoretical works have reported how the mismatch 
strain between the precipitate and the surrounding ma- 
trix phase changes the characteristics of the microstruc- 

tureiEiiaiiaaiaiiE3 

One of the primary reasons for considering the role of 
elasticity in microstructure formation is that it offers us 
the possibility of using elastic stresses to design desirable 
microstructures. In this context, one of the possibilities 
often cited is that of inverse coarsening. In inverse coars- 
ening, the elastic fields promote the growth of a small 
precipitate at the expense of a larger neighbor - exactly 
the opposite to surface energy driven coarsening. In prin- 
ciple, elasticity could favor an array of equally sized pre- 
cipitates, which would lead to the improved mechanical 
properties of the alloy. 

Particle splitting is perhaps the most closely linked 
elastic effect to inverse coarsening. Particle splitting is 
a process in which one cuboidal particle splits into sev- 
eral particles, usually two (a doublet) or eight (an octet). 
Such splitting has been observed experimentally 0, y| 
and has been theoretically identified as an elastically 
driven process 0, . 

While splitting may be an important factor for inverse 
coarsening, its mechanism and the precise experimental 
and theoretical conditions under which it occurs are not 
completely clear. Some experimental pictures suggest 
that splitting occurs via a morphological instability in 
which the sides of a cuboidal precipitate become increas- 
ingly concave prior to splitting [3. Other pictures sug- 
gest a hollowing at the center of the particle [H. Many 
theoretical simulations of splitting have used diffuse in- 
terface techniques, and suggest that the instability oc- 
curs by the 'hollowing' mechanism However, the 



hollowing mechanism cannot explain the formation of an 
octet and the fact that the splitting occurs mostly in 
the concave particle shape 0, 1^ though it can describe 
splitting to a doublet. Recently, Lee, using a discrete 
atom method, observed splitting by an instability mech- 
anism IqI. 



In this letter we demonstrate, through phase field 
model calculations, that interface instability driven by 
elastic anisotropy and diffusion field can generate parti- 
cle splitting. Splitting occurred only under the concave 
growth condition and the picture that a concentrated 
elastic field around interface groove promotes splitting 
is very similar to Asaro-Tiller-Grinfeld (ATG) instabil- 
ity 1 ij . In addition, this splitting mechanism can explain 
the formation of both doublets and octets. 

The model in this study assumes coherent interfaces 
(i.e. there are no dislocations) and takes into account 
the effect of anisotropic elasticity and elastic inhomo- 
geneity which represents that the elastic constants of the 
matrix and the precipitate differ. We assume that the 
misfit between the matrix and the precipitate is purely 
dilatational with magnitude e^, the elastic constants of 
both phases have cubic (four-fold) symmetry, and that 
the lattice parameter is not a function of composition. 
The model is formulated in terms of the concentration 
field c{f,t), the phase field (p{r,t) which is zero in the 
matrix phase and one in the precipitate, and the elas- 
tic strain tensor — 1/2 {dui/dxj + duj /dxi), where 
M(f, t) is the atomic displacement field. The interfacial 
region is defined to be a mixture of matrix and precipitate 
with a different composition but with the same chemical 
potential Therefore, the concentration field of the 

system is defined as c = c^h{(f)) + c'^^[l — /i((/i)], where 
is a concentration in the A-phase and h{(j)) is the smooth 
interpolation function, which is zero &i (p — Q and one 
at = 1 and changes smoothly from zero to one, i.e. 
< (/) < 1. The superscript P and M represent the pre- 
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cipitate and matrix phases, respectively, which will be 
used later without description. From the definition of 
the interfacial region, the following constraint should be 
satisfied at any point in the system: 



(1) 



where / is the chemical free energy of the A-phase and 
represents the derivative of with respect to c^. 
The free energy functional is: 



df 



(2) 



where /(</), c, Uij) is the free energy density of the system 
including the coherent strain energy, defined as 'wg{4i) + 
/^(c^)/i(0) + /^(c^)[l - + /^'(</., Uij), where g{^) 
is the double well potential, which has minima at 1 and 
0, the equilibrium values of the precipitates and matrix 
phases, respectively. In this study, g{(j)) was selected as 
— (j))^. f^^ is the coherent elastic energy density 
caused by the lattice mismatch between two phases. In 
general, the elastic energy density is expressed as [l^ : 



r 



-Ci-jkii(f>)[uij - e^d^jip{(f))][uki - e'^Skitl'i(t>)] (3) 



where Cijki is the stiffness tensor that depends on the 
phase field due to elastic inhomogeneity and ipi^t') is a- 
smooth function with the same characteristics as h{(j)). A 
summation convention over repeated indices is implicit. 
Considering the spatial inhomogeneity of the elastic stiff- 
ness, the stiffness tensor is given by: 



(4) 



where ^{(f>) is a smooth function which is zero at = 1 
and one at </> = 0, Cf^,^ = C^,, , and = C^,^ - 
where C^^i is the elastic stiffness of the A-phase. 

It is reasonable to suppose that the elastic field relaxes 
much faster than c or 0. The elastic field can then be 
solved in terms of the phase field using the condition of 
local mechanical equilibrium: 



oui 







where aij = SJ- /Suij is the stress tensor: 

<Jij = Cyfci(0)[ufci - e^4;V(0)] 



(5) 



(6) 



In this work, in order to obtain accurate solutions of 
equation Q with minimal calculation time, the high or- 
der approximation method similar to the method pro- 
posed by Chen and coworkers is used instead of the 
iteration method which gives very accurate solutions but 
requires considerable calculation time. 

The solution of equation © to zeroth order in the 
elastic stiffness C^'^-^,; is 



(7) 



where A represents the value of A in Fourier space, qj 
the j-direction component of the wave number vector. 



e^Skh and a^/, indicates a component of the 
f aik which is defined as Cfj^iljl^ 
the displacement field to the nth order becomes 



inverse tensor of aik which is defined as Cf^kiljl^- Then, 



^k + i(^kl(^i3lme^^lynq]{i{4')i'{4')}l 



(8) 



where vJl~^ is the displacement field to the n-ith order 
and {}f represent the Fourier transform of the term in 
braces. 

The elastic field can now be expressed in terms of the 
phase field 0. Substituting the solution for the strain field 
gives the free energy in terms of c and 0. Assuming that 
the system evolves in time so that its total free energy 
decreases monotonically, the evolution equation for the 
phase field becomes: 



dt 



0(p 



(9) 



with 



U = h'mr-f'' + {c''-nf'M + ff+wg'{(^) (10) 
while the concentration field should evolve with: 



dt 



= V 



fee 5 c 



= V ■ D{(t))Vc + V 



Ice 



fe^y4> (11) 



where M is the phase- field mobility, £>((/)) is the diffusiv- 
ity of solute atoms, and f^, fee, and fc^ are df/d4>, 
df^^ /d(f), d"^ f ld(? , and d'^f/dcdcj), respectively. 

In the sharp interface limit, the phase field model pre- 
sented here reduces to a modified Gibbs-Thompson equa- 
tion with a coherent interface as follows: 



aavn — KeCr + A/ip + G'^ 



(12) 



Here, a — D/Me^, Vn is the interface velocity, a the in- 
terface energy, Kc the curvature of the interface, and A/ic 
the chemical driving force of the phase transition defined 
as fiQ — /ip^, where represents the chemical potential 
at the A-phase side of the interface. 6"^°'* means the elas- 
tic energy required to maintain coherence at the interface, 
defined by (uf- — ufAaff , where uf, and (t,4 are elastic 

•J \ LJ ''J ' ''J f-J ''J 

strain and stress tensors at the A-phase side of the inter- 
face, respectively 0]. G*'"""'' represents the difference in 
the elastic strain energies of the precipitate and the ma- 
trix, as defined by l/2a,^(u^ - e^Jy ) - l/2affuff 
The detailed derivation of Eq. (|12l) will emerge else- 
where [l^ . 

Numerical simulations on a discrete lattice were per- 
formed in two-dimensions. Euler's method was used for 
the integration in time. For all simulations presented 
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here, the mesh size Aa; = 1.0 x 10 ^m, the time step 
At = 0.0125, and e = 3.52 x 10"^ This choice of Ax 
and e guarantees that the interface is resolved by at least 
five points and the interfacial energy is 0.0243 J/m? . We 
assume that the diffusivities of two phases are the same 
and 1 X lO^^^m^/s. Quadratic functions were used for the 
free energy functions of both phases for simphcity, that is, 
fP = A^(c^ - clflVm and = A^(c^ - cf )VK., 
where A^ and A^^ were 10000 and 2000, respectively, Vra 
is molar volume and was set at 7.0 x 10~^to^/too^, and 
and are the equilibrium concentrations without elas- 
tic field and were set at 0.2 and 0.4, respectively. The 
elastic constants of two phases (in units of VS^^ er g j crn?^ 
are Cfi = 1.67, Cfa = 1.07, Cj^ = 0.99, = 1.12, 
Cli = 0.63, and Cf^ = 0.57, respectively; for the Ni- 
based superalloy j^. The misfit strain between the pre- 
cipitate and the matrix is dilatational and assumed to 
be 0.8%. The displacement field was calculated to the 
third order from Eq. (jHJ because the inhomogeneities of 
the elastic constants were very large. Periodic boundary 
conditions were employed in all directions. 
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FIG. 1: Sequential evolutions of the precipitate (the top row) 
and the corresponding misfit elastic energy density distribu- 
tions(the second row): from the left, t = 3.125, 31.25, 62.5 and 
150s. The lowest row shows splitting and split 7/-precipitates 
in the Ni-based superalloy system Q. In the calculation, 
M = 5.0 x 10"^ The total size of the calculation domain 
is 2.56fim x 2.56/^m. 

Figure n shows the calculated splitting behavior of the 
precipitate and the corresponding elastic energy density 
distribution at various times. The morphological in- 
stability of the interface due to the diffusion field and 
anisotropy of the elastic stiffness is evident and this in- 
terface instability leads to particle splitting. Despite the 
simplicity of the two-dimensional calculation, the split- 
ting pattern is in good agreement with the experimental 



results (see Fig. and j). The scenario of the splitting 
is as follows: initially the precipitate grows in a circu- 
lar shape due to the isotropic interface energy, and then 
transforms to cuboidal shape above a critical size, be- 
cause the misfit strain energy proportional to its vol- 
ume dominates the interfacial energy in proportion to 
the interfacial area. Subsequently, the growth rate in 
the < 11 > direction is enhanced due to a high chemi- 
cal driving force while the migration rate in the < 10 > 
direction is suppressed due to enriched solute and a cor- 
responding low chemical driving force. Therefore, the 
growing shape transforms to a concave shape. In the 
concave growth regime, the interface in the < 10 > di- 
rection has a negative curvature, and so, the elastic strain 
energy concentrates around the concave edge . This con- 
centrated elastic energy drives the dissolution of the pre- 
cipitate, and then a groove forms and penetrates the pre- 
cipitate (see Fig. ^). This situation resembles Asaro- 
Tiller-Grinfcld instability, in which elastic media under 
the nonhydrostatic stress show the morphological insta- 
bility [11|. Through various numerical simulations, we 
confirmed that the splitting occurs only in the concave 
growth condition and that it is induced by the interface 
instability. This is in good agreement with the experi- 
mental observations 1^ • 
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FIG. 2: Time variations of < 10 > length (ai) and < 11 > 
length (a2) of the precipitate (a) and the corresponding ratio 
of ai to a2 (b) which is defined as the shape parameter C,. 
The inset shows the definitions of ai and 02. In the case of 
splitting, M = 5.0 X 10^^ and in non-splitting, M = 3.0 x 
10-^ 

Figure 121 shows the sequential variations of the < 10 > 
length (ai) and < 11 > length (02) of the precipitate and 
the ratio of the two lengths which is called a shape param- 
eter (C) hereafter. There are four different growth modes 
(see Fig.|21b). In the first growth mode the growth rates 
of ai and as arc proportional to i0.35±o.oo5 ^^^^ ^o.55±o.oi^ 

respectively and hence transition from circular through 
cuboidal to concave shapes occurs (for a perfect square 
C — 0.7). The second regime can be called the shape- 
conserving mode because the shape parameter remains 
constant in this period. The growth rates of ai and 02 
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have the same dependence i0-3i±0-0i on time. The growth 
of oi stops upon entering the third growth mode while 
02 continues to grow and the shape parameter starts to 
decrease again. In this regime the grooves are gener- 
ated on the four interfaces of the < 10 > direction (see 
Fig. and concavity around the grooves starts to relax 
because the elastic field around the grooves is relieved 
by local concentration of the stress field on the edges 
of the grooves. One can see this by comparing Fig. 
with Fig. m. Once the elastic energy concentrated on 
the tips dominates the chemical driving force for precipi- 
tate growth which is reduced due to the buildup of solute 
around the grooves, the grooves penetrate into the pre- 
cipitate and the fourth growth stage begins. Note that 
the particle in the non-splitting case also shows concave 
growth, but then returns to the cuboidal shape with time 
(see the variations of the filled triangles in Fig.[2l3). 



ered as a meta-stable phase and the activation barrier 
is the interfacial free energy corresponding to the in- 
creased interfacial area. Considering these points, the 
splitting phenomenon is kinetically rather than energet- 
ically driven. 

In this work, it has been shown through phase field 
study that the interaction between the elastic field con- 
densed at the concave region and solute enrichment 
causes the particle splitting phenomenon. The splitting is 
mediated by interface instability, which differs from pre- 
vious studies that splitting occurs through the nucleation 
of the matrix phase at the center of the precipita te M, 
or that it is caused by energy considerations Il0l|. In 
particular, the sequential splitting behavior realized by 
our simulation shows an excellent agreement with previ- 
ous experimental studies 0, 1^ • 
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FIG. 3: Time variations of the total free energy (Eq.|2j, elastic 
energy (a), and the total free energy, except the elastic energy 
(b). The phase field mobilities are 5.0 x 10~^ for the splitting 
case and 3.0 x lO"'^ for the non-splitting case, respectively. 

Fig.Oshows the variations of the total free energy [Eq. 
(PI], elastic energy, and the total free energy excluding 
the elastic energy, with the evolution of the second par- 
ticle. Here, the third free energy term mentioned above 
is called the chemical free energy although it includes 
the interfacial energy (i.e. gradient energy and double 
well potential energy). During the evolution of the par- 
ticle, the total free energy and the chemical free energy 
decrease while the strain energy increases due to misfit- 
ting particle growth. When the splitting has been al- 
most completed, the chemical free energy jumps due to 
the increased interfacial area, but the total strain energy 
reduces by interface grooving, which actually resembles 
the ATG instability [ll|. As the decrease of strain en- 
ergy overwhelms the increased interfacial free energy, the 
total free energy decreases and the splitting process ad- 
vances. As shown in Fig. the configuration of four 
split pieces is a more stable state than non-split single 
particle. Hence the non-split large particle is consid- 
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